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We study the Kohn-Luttinger effect in a two-dimensional (2D) nested Fermion iiquid with a re- 
pulsive interaction via the renormalization group method and identify the resulting order parameter 
symmetry. Using the band structure of the 2D Hubbard model close to half-filling as a prototype, 
we construct an effective low-energy theory. We use multidimensional bosonization to incorporate 
the zero sound channel and find marginal Fermi liquid behavior in the absence of any instability. We 
show an analog of the Landau theorem in nested Fermion liquids, which serves as the criterion of the 
BCS instability. Including repulsive or antiferromagnetic exchange interactions in the low-energy 
theory, we show that the d x 2 _ y 2-wave BCS channel is renormalized to be the most attractive. Be- 
low half-filling, when the nesting is not perfect, there is competition between the spin-density-wave 
(SDW) and the BCS channels; when the SDW coupling is small enough, there occurs a d a ,2_ H 2-wave 
superconducting instability at sufficiently low temperatures. 
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I. INTRODUCTION 

One question of recent interest concerns the possible breakdown of the Fermi liquid in the presence of singular or 
even regular interactions. One example is the Luttinger liquid in one spatial dimension where the quasiparticle pole is 
destroyed by a short-range interaction although in this case the low-energy excitation remains gapless. The existence 
of the Luttinger liquid in higher spatial dimensions has also been proposed but is still controversial .El Another example 
is the Kohn-Luttinger (KL) effect: A Fermi liquid state is unstable against one or more BCS channels in the presence 
of a repulsive interaction. □ It has been shown that a short-range repulsive interaction can induce an attractive BCS 
channel at large orbital angular momentum modes. In three-dimensional Fermi gases, the Kohn-Luttinger effect was 
found in the p-wave BCS channel^ This notion was extended to two-dimensional Fermi liquids and a weaker KL effect 
was also found in the p-wave channel although the effect is not visible in a second-order perturbative expansion. □ 

However, the KL effects in isotropic Fermi liquids are weak because they are found in sub-leading corrections in the 
renormalization group (RG) sense; the resulting T c is expected to be low. For example, in three dimensions, the KL 
correction to the function is found to.be 0(s~ 7 ^ 4 ) where s is the dimensionless length scale which is taken to be 
a large number in the low-energy limit B But in the case of nested Fermion liquids, pthe effect is more dramatic since 
the nesting gives rise to a logarithmic flow in the repulsive channel. Baranov et aln studied the KL effect in a two- 
dimensional Hubbard model at low filling factors where the Fermi surface is almost circular and found an attractive 
channel with d^-wave symmetry, but it is also a sub-leading effect which originates from the directional property of 
the square lattice. In this paper we study the KL effect in two dimensions in a nested square Fermi surface. 

The subject of this paper may be of interest in the context of the cuprate high-T c superconducting materials. 
The KL effect implies that an initially repulsive interaction between electrons can induce an effective attractive BCS 
interaction among quasiparticles lying in the thin momentum shell aroruad the Fermi surface. Hence this is one possible 
candidate for a non-phononic mechanism of cuprate superconductivity.!] Indeed, we find that the repulsive interaction, 
which naturally contains an antiferromagnetic (AF) exchange channel, favors a d x 2_ y 2-w&ve BCS instability. This is 
in accordance with expectations Eel 

The low-energy theory of the two-dimensional Hubbard model is complicated by several features of the anisotropic 
nested Fermi surface. Especially near half-filling, there are van Hove singularities in addition to non-uniformity of 
the Fermi velocities along the Fermi surface. Moreover, in the low-energy limit, the system is not a regular Fermi 
liquid due to nesting. For instance, the quasiparticles are expected to have a short lifetime, 1/r ~ T, near the Fermi 
surface, and the existence of a well-defined Fermi surface may be questioned.S To incorporate these characteristics, 
we introduce a low-energy theory which models the bancLsteucture of the two-dimensional Hubbard model near 
half-filling. We make use of multidimensional bosonizationHi!3 to describe the low-energy properties of the model. 
Although Fermi liquid theory may break down near the Fermi surface, marginal deviations from a Fermi liquid are 
within the scope of bosonization. 

The RG in the bosonic basis is a suitable method to study the low-energy phases of the model since it puts all 
interaction channels on an equal footing. Wilson's RG algorithm, which we shall make use of, consists of integrating 
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out high-energy degrees of freedom to construct a low-energy effective theory. In a degenerate fermion system, this 
amounts to successively eLkpi&ating inner and outer sides of the momentum shell around the Fermi surface in which 
the low-energy theory liesBll^llJ In this paper we construct a low-energy effective theory of a nested fermion liquid by 
the RG procedure and then investigate the RG flow of the effective interactions in the framework of bosonization. In 
other words, we successively integrate out higher-energy modes of the bosonic field <p. In most cases, this bosonized 
RG procedure yields results which agree with those obtained in the fermion basis. The distinct advantage of the 
RG treatment based on bosonization is that the existence of well-defined quasiparticles is not a priori assumed. 
Furthermore, those interaction channels which preserve the infinite-/7(l) symmetry^ on the Fermi surface are exactly 
diagonalized at the outset £j 

This paper is organized as follows. In Sec. [n], we introduce the model low-energy theory of a two-dimensional 
nested Fermion liquid. We use multidimensional bosonization to incorporate the zero sound channel and show that 



the resulting low-energy fixed point is a marginal Fermi liquid. In Sec. [II, we derive the RG equations for the 
BCS channel and discuss the condition of the BCS instability. An analog of the Landau theorem is applicable to 
the case of anisotropic Fermi surfaces. In Sec. IV, we incorporate marginally relevant repulsive or AF interactions 
and evaluate their RG flow which contribute to the BCS channel. The d x 2_ y 2-wave channel is the most important 
one, and we argue that below half-filling competition arises between the spin-density-wave (SDW) and d a ,2_ y 2-wave 
superconductivity. A brief summary and conclusions are given in Sec. |v|. 



II. MODEL LOW-ENERGY THEORY 



The Hubbard model is a convenient description of strongly correlated fermion systems. For a reliable investigation 
into the RG of this model, it is essential to accommodate a properly realistic band structure near the Fermi surface. 
Here we present a model low-energy theory which simulates the nested fermion system. 

First, consider the band structure of the Hubbard model in two dimensions near half-filling. The Hamiltonian is 

H = -t^24 a c ja + U22mi mi , (I) 

where i,j are between nearest neighbors on a square lattice and U is the on-site interaction strength. Here c\ a and 
Ci a are fermion creation and annihilation operators. Throughout this paper we assume that \U/t\ <C I holds so that 
the perturbative expansions are justified. For now we turn off the on-site interaction term and discuss the character 
of the Fermi surface that results from the hopping term. The spectrum of the fermions is 

E(k) — — 2i(cos k x a + cos k y a), (2) 

where a is the lattice spacing. At half-filling, the Fermi surface is described by the points in k space which satisfy 

cos k x a + cos k y a = . (3) 

The Fermi velocity is a function of k given by 

v(k) = 2ta(x sin(fc x a) + y sin(fc y a)) (4) 

on the Fermi surface. Singularities in the density of states occur at four corners of the Fermi surface, k = (±ir/a, 0) 
or (0, ±7r/a), where the Fermi velocity vanishes. n 

Turning on the on-site interaction, the RG may be used to integrate out the high-energy degrees of freedoms The 
resulting low-energy effective theory is expressed in terms of quasiparticles living in a narrow shell around the Fermi 
surface of energy cutoff s c . Quasiparticles may then be related to the bare fermion operators by 

V^k = Zy. 1/2 (e c ) c k , (5) 

where the wave-function renormalization factor Zk(E c ) rescales the discontinuity in the quasiparticle occupancy at the 
Fermi surface back to 1. The final quasiparticle weight Zk can be obtained by further integrating out the remaining 
low-energy degrees of freedom, lim £c _ + o Z^(e c ) = Z^. In non-Fermi liquids, Z^ is zero. Here we stop the RG process 
at the energy scale e c and investigate the low-energy theory thus obtained. Assuming that the weak nearest-neighbor 
interaction term does not alter the shape of the Fermi surface or the band structure in Eq. (||) , we may linearize the 
spectrum of the quasiparticles: 
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E(k) =v \sm[k F (k) x a}\ g y (k) , (6) 

which is similar to Eq. (||) close to the Fermi surface. Here qn (k) = (k — kp(k)) • n(k), where kp(k) is the wave- vector 
on the Fermi surface closest to k, k F (k) x is the x component of kp(k), and n(k) is the unit vector normal to the 
Fermi surface at kp(k). 

The difficulty with bosonizing this Fermi surface is that at the four corners of the Fermi surface the Fermi velocities 
are zero and the Fermi surface normal vectors are not well-defined in the vicinity of the van Hove points. Below 
half-filling, however, the van Hove singularity is not encountered. So we introduce a model spectrum which simulates 
the Hubbard model near half-filling but with a smoothed out van Hove singularity. We modify the spectrum in Eq. 
(§) as 

E(k) = v \sm[k F (k) x a}\ gy(k) , e/a < \k F (k) x \ < (tt - e)/a , 

= v sin(ea) qn (k) (otherwise) (7) 

and we smooth the four corners of the Fermi surface as four quarter-circles with a radius \/2e/a (Fig. [l]). So we 
define the Fermi velocity function v F (k) as the coefficient of <?||(k) in Eq. (fjj). This choice is reasonable because the 
number of states near the van Hove points here is of order e/sin(e) ~ 0(1), which is finite even if the limit e — > 
is taken. This artificial smoothing simulates the van Hove singularity as the van Hove points accommodate a fixed 
fraction of the total number of states. From now on we set 1/a = 1 and v — 1 for convenience. The factor 1/a may 
be recovered whenever a comparison with other momentum scales is needed. With this convention, for example, the 
nesting vectors are expressed as Q = ±(n, tt) or ±(n, — n). 

Now that the Fermi velocities are well-defined over the whole Fermi surface, the model may be bosonized. We coarse- 
grain the Fermi surface into small squat boxes of a width A along the Fermi surface and a height A(S) = £ c /vf{S) 
normal to the Fermi surface, as shown in Fig. |l|. Here v F (S) is the Fermi velocity function at the squat box S. Then 
we introduce a coarse-grained charge current which lives in a squat box labeled as S: 

J a (k s ; q) = 0(S; k - q) 9(S; k) {^_ q Q Vk,« - 5^ n k , a } . (8) 

k 

Here rik,a = a ipk,a)- & is the spin index and the label S identifies a patch on the Fermi surface with a Fermi 
momentum ks and 9(S; k) = 1 if k lies inside a squat box centered on ks and 9(S; k) = otherwise. The two 
momentum scales A, X(S) are taken to be small so that 1/a » A » ■MS) & t eacn Fermi surface patch. The latter 
inequality is needed to minimize scattering of fermions outside of the squat box so that the number of fermions in 
each patch on the Fermi surface is conserved. Now the limit of e — > is not permitted since it violates the condition 
of e/a 3> A at the four van Hove opiats. Instead, we take e as a small non-zero number. 
The U(l) current algebra readsETtil 

[J a (k s ; q), Jp{k T ; p)] = ^^S.T^q+p.o^ q ' n ~s , (9) 

which holds when 1/a > A » A(S). Here fl = A(L/2tt) 2 is the number of states in the squat box divided by A(S) 
and rig is a unit vector normal to the Fermi surface at patch S. The algebra, Eq. (^), makes it possible to bosonize 
the currents. The bosonic representation of the currents is 

J Q (k s ;x) = V4^ n s • V0 a (k s ;x) , (10) 

where the bosonic fields <j) a satisfy the canonical commutation relations. Then the fermion quasiparticlc fields at each 
patch can be expressed in terms of boson fields asE£l 

^ Q (ks; x, t) = e* ks ' x J2 *(S; pK (p - ks) ' x Vw (*) 



l A KS)/2 iks . r.-VSF 

(2tt)2 6XP 



{^^(k s ;x,f)}6(S). (11) 



Here O(S) is an ordering operator introduced to maintain Fermi statistics along the Fermi surface. 
The bosonized effective action of the low-energy theory is 

S = J2 E [ ^-v F (S)q-h s ]a* a (k s ;q) a a (k s ;q) + S I . (12) 

S q,A(S)>qn s >0 ^ 
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Here Si is the fermion interaction terms which we discuss later. The charge currents are related to the canonical 
boson operators a and a' by 



J Q (k s ;q) = \/£l\q- n s | [a a (k s ;q) 6»(q • n s ) + a^(k s ; -q) 



-q - n s J 



(13) 



where 



[a a (k s ;q), at(k T ; p)] = SapSs^S, 



2 

p,q 



(14) 



Here 0{x) = 1 if x > and 0(x) = otherwise. From Eq. ( fLl| ) the fermion. Green's function can be expressed in terms 
of the Fourier-transformed boson correlation function and it is given bylij 



G£(k s ;x,<) = i (V tQ (k s ; x, t) ^ a (k s ; 0, 0)> 



x j^l (a a (k s ;q) oi(ks;q)) 
Aq ■ n s 



^.r p i(q-x-wt) 

2tt L 



(15) 



Next consider the renormalization of the fermion interactions. As in the case of a circular Fermi surface, the marginal 
interactions take the form of four- fermion (two-body) interactions, as can be verified by counting the engineering 
dimension. Unlike the case of a circular Fermi surface, however, many large-momentum scattering channels are now 
marginal due to the geometry of the nesting. Some of the channels give rise to spin-density-wave (SDW) or charge- 
density- wave (CDW) instabilities. These logarithmically renormalized channels are discussed in Sec. [V. For now, 
the zero-sound (ZS) channel may be incorporated. This channel yields a Fermi liquid solution for a circular Fermi 
surface in two dimensions and a Luttingcr liquid in one dimension. It is important to include this channel in the 
free action before the the large-momentum scattering channels are studied perturbatively. This permits a check on 
whether the Fermi surface breaks down in the presence of the ZS channel. 

For simplicity, we suppress all momentum dependence in the coupling function and drop the spin indices. Inclusion 
of the spin exchange channel does not change the result qualitatively. The ZS channel in the fermion basis is 



Szs — 



d 2 x dt F(S,T) ^(ks;x,t) V(ks;x,t) ^(k T ;x,t) V(k T ;x,t) . 



(16) 



Note that a U(l) phase rotation in each squat box, t/>(kg;x, t) — > e ir ^ s '-0(ks; x, t) and ^t(kg;x, t) — > 
e~ ir ^ip^ (ks; x, t), leaves the action invariant. Physically, this infinite U(l) symmetry means that the number of 
fermions at each patch on the Fermi surface is conserved. This channel can be expressed in terms of U(l) currents 
since Eq. (p|) is U(l) invariant at each patch. The ZS channel in the bosonized form is succinctly written as 



~>zs 



S,T q 



2^ 



F(S,T) J(k s ;g) J(k T ;-g) , 



(17) 



which is marginal since it is bilinear in bosons. Here q denotes (q, to). The solution is a Gaussian boson parametrized 
by the shape of the Fermi surface and the coupling parameter F(S,T). Those which do not have infinite U(l) 
symmetry are not expressed as bilinear in currents but they can be studied perturbatively. Here for simplicity the 
case F(S, T) = F is studied. 

Turning our attention to the fermions on the nested part of the Ferrai .surface, since that is where the most singular 
behavior may arise, we make use of the generating functional mcthod,u3ii3 and obtain the boson correlation function 



(a(k s ;q)a t (k s ;g)) = 



A 



q n s 



~ v F (S)d-& s (2tt) 2 (u-v F (S)q-n s )2 
and the fermion Green's function is given by Eq. (|l5|). Here, 



(18) 



K (q) 



1 
1 



2A 

(27T) 



E 



^(q-n s ) q- n s 

lo - v P (S)q ■ n s 



n r 

— \{2tt) 2 xo(^,sine q) +XN(w,q) 



(19) 
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where xn is the contribution from the nested part of the Fermi surface: 



xn{q) 



In 



1 — uJ 2 /q 2 + cose- 



/gji — cose 



■An 



i/l — uj 2 /qj_ + cos e 



\/l — u 2 /qj_ — cos e 



(20) 



Again q» is the component of q parallel with the Fermi surface normal and q± is the perpendicular component. The 
contribution from the curved-part of the Fermi surface (xo) m Eq. ( P^l) gives a self-energy correction Im S ~ u 2 In ui 
as in a Fermi liquid system.Ej We perturbatively estimate the imaginary part of the quasiparticle self-energy which 
originates from Eq. (|20|). The imaginary part of xat(q) arises when \lu\ < \q^\ sine or \w\ < \q±\ sine. To second 
order in F , the boson correlation function is expanded as 



i(k s ;g)a t (k s ;g)) 



+i 



v F (S)q • n s 

A q • n s 



(2tt) 2 (u-v F (S)q-n S ) 



F 



\/2 

F 7^lY2 XN(q) 



(2tt) 2 



(21) 



where the xo term is dropped. The first order term in Fq renormalizes the Fermi velocity at the quasiparticle pole 
just as in the Fermi liquid case, but the velocity correction is dv = FoA/kp(2ir) 2 , which is infinitesimal. Spin-charge 
velocity separation does not occur, as can be easily shown when the spin indices and the exchange channel are 
included in the effective action. Inserting Eq. ( pi] ) into Eq. (15), expanding the exponential to second order in 
and performing a Fourier transform into (k, u>) space yields the imaginary copiection to the quasiparticle self-energy.EEl 
Then Im S(k, ui) cx lu F^k/kF, which is of a marginal Fermi liquid form.Ea Similarly, corrections of higher orders 
in Fq can be estimated and we find that they also yield a marginal Fermi liquid form. This agrees with previous 
expectations by others who obtained the inverse lifetime of quasiparticles, 1/r ~ T, as a function of the temperature.tJ 
Since the quasiparticle weight of a marginal Fermi liquid vanishes at the Fermi surface as the inverse of a logarithm, 
the Fermi liquid fixed point breaks down but the breakdown of the Fermi surface is smooth enough for the bosonization 
and the subsequent RG procedure to be applicable £3 Therefore we can confidently apply the RG method to analyze 
large-momentum channels which do not preserve infinite {/(l)-symmetry. 

One might wonder whether the marginal- Fermi- liquid- like behavior of the., nested fermion system has relevance 
to the linear-T in-plane resistivity of cuprate materials in the normal state.113 However, we only have considered a 
simplified single-band Hubbard model with the ZS channel alone whereas a three-band Hubbard model may be a 
better description of the cuprate materials. For example, .-we have not incorporated in the model the strong van 
Hove singularity which is present in the cuprate materials.c3 Furthermore, the n-doped N2_ a; Ce 2: CuO 4-5 (NCCO) 
material exhibits T 2 resistivity even though photoemission studies show a nested Fermi surface.Eil This suggests that 
nesting with the ZS alone is not sufficient to account for the behavior of the resistivity in the normal states of cuprate 
superconductors. 

In the following sections, we study large-momentum scattering channels by the perturbative RG method. We do 
not explicitly include the ZS channel at the outset but we assume that it contributes sub-leading corrections to RG 
equations. For example, though bosons in different patches are correlated via the ZS channel, the correlation between 
fermions in different patches is negligible: 



(V> T (k s ;x) ^(k T ;0)) cx -e 



<0(k s ;x)0(k T ;O)) 







(22) 



where L is the system sizeJli Also the self-energy corrections are of order Act <C 1 times smaller than the leading 
order. So the leading RG equations we obtain in the following sections are assumed to be unchanged even if the ZS 
channel is considered. 



III. RG OF THE BCS CHANNEL 



In this section we study the RG flow of the BCS channel for the nested Fermi surface and find an analog of the 
Landau theorem. As for the case of an isotropic circular Fermi surface, the BCS channel in a nested square Fermi 
surface flows logarithmically in scale length since the Fermi surface is symmetric under reflection at the origin in 
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momentum space. On a circular Fermi surface, the Landau theoremEl states that there occurs a BCS instability if 
any one value of Vi is negative where Vi is the Ith orbital angular momentum component of the BCS vertex function 
V{di,02) — Yli e ll ( 9l ~ e ^ Vi- Here 01,6*2 are the angular variables aa-.the Fermi surface in momentum space. This 
theorem can be seen in the RG equation for Vi derived by Shankarjjli3 



the solution of which is 



dVi 
d In; 



Vi(s) 



2vr 1 



V, 



(23) 



(24) 



for each I. Thus, if any Vi < 0, it eventually leads to instability. 

In the case of the nested Fermi surface thera-is no rotational invariance in the BCS vertex function, and so T^(ks, kx) 
must be expanded in a double Fourier series .3 Consequently the Fourier coefficients are coupled in the resulting RG 
equations, unlike in the case of a circular Fermi surface. First we discuss the possible symmetry of the order parameter 
in a square lattice. The BCS vertex function has D4 symmetry and there are five sets of basis functions which generate 
the two-dimensional function space to which the. vertex function belongs. The five sets are represented as A\ and A 2 
(s-wave), B\ and B 2 (d-wave), and E (p-wave).0 The basis functions in each set are following: 
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B 2 
E 



cos 2nk x , 

sgn(fcj,) sin2nk x , 

cos(2n + l)k x , 

sgn(fcy) sm(2n+l)k x , 

a sin(n + l/2)k x + b sgn(k y ) cos(n 



l/2)k x 



(25) 



where we have parametrized the basis functions with k x and the sign of k y on the square Fermi surface which is defined 
by cos k x + cos k y = 0. In general the BCS vertex function on a square lattice can be expressed as a double Fourier 
expansion in terms of these basis functions. Due to the D4 symmetry of the BCS vertex function, the expansion does 
not contain crossed terms between basis functions of different symmetry sets, and so the five sets are decoupled from 
one another in the RG equation at the outset. !_. 
Now we inspect the RG equation of the BCS obtained in bosonization language,E3 

^ (kS ' kT) - A E^^-M^.M. (26) 



d(lns) 



(27T) 



Notice that there is a weight function of 1/vf(U) in the sum. To analyze the RG equation, we modify the basis 
functions in Eq. (|2^) so that they are orthonormal under a summation over the Fermi surface with the weight 
function l/vp{U). The modified basis functions are given by 



ai(n,k) 



a 2 (n,k) 



&i(n,k) 



& 2 (n,k) 



2ttV2 
sgn(ky) 

_2nV2 
sgn(ky) 
t'F(k) 



1/2 



cos2nk x , ai(0,k) 



wf(k) 



4vrV2 



1/2 



1/2 



2ttV2 

1/2 

cos(2n 



sin 2nkr, 



l)k x , 



2ttV2 



2ttV2 

1/2 

h 



1/2 



e(n, k) 

The BCS vertex function can be expanded as 

V(k s ,-k T ) = J2{c? m oi(n,ks) ai(m,k T ) 



sin(2n + l)k x 



sin(n + l/2)k x + C2 sgn(fc;,) cos(n + l/2)fe x ] 



(27) 



C 2 "™ a 2 (n,k s ) a 2 (m,k T ) 



+Cr i &i(n,k s ) &x(m,k T ) + CT 6 2 (n,k s ) b 2 (m,k T ) 
+C£ m e(n,ks) e(m,k T )} . 



(28) 



G 



This expansion is valid as long as minf|uir(k)|^ > which holds in the model considered here. In fact, min^|wi?(k)|^ 

is zero only at exactly half-filling. 

Now insert Eq. (^8|) into Eq. (|26| ) to obtain the RG equations of the Fourier coefficients: 

dE7 = -c^£c?'3 m (29) 

v ; 1=0 

where i = 1, 5. Now the analog of the Landau theorem in the nested fermion liquid can be obtained. Namely, the 
diagonal Fourier coefficients obey the equations 

§7=-(^I>;r. (30) 

v ; 1=0 

Here the diagonal coefficients decrease as the length scale increases. Thus, if any one value of C™" starts out negative, 
a BCS instability occurs. This is similar to the Landau theorem in an isotropic circular Fermi surface. This criterion 
of instability permits us to identify the possible symmetry of the resulting order parameter, since the gap function is 
related to the BCS vertex function by the following: 

A k = -Enk,k ^, (3D 

where k, k' are within the thin shell around the Fermi surface and E(k!) is the spectrum of the low-lying excitation 
in the superconducting state. In the next section we apply this result to search for a possible BCS instability in the 
presence of a repulsive interaction. 



IV. RG OF REPULSIVE INTERACTIONS 



In this section we discuss the RG flow of repulsive interactions. We only consider those channels which flow 
logarithmically in length scale s. To simplify the analysis we take the coupling functions to be independent of 
momentum. 

There are three classes of interaction channels which are logarithmically renormalized in the presence of nesting. 
The first class is the CDW or SDW channel. In case of a repulsive bare Hubbard interaction, it is the SDW channel 
that leads to an instability. This channel is shown to have little overlap with the BCS channel and is practically 
decoupled from the BCS channel. The second and third classes are Q = and Q = 2fc^ repulsive (or AF exchange) 
channels which renormalize the BCS channel. In the following two subsections we study the RG equations in the 
SDW (CDW) channel and in the other repulsive (or AF) channel which renormalize the BCS and calculate the 
Kohn-Luttinger effect. 



A. CDW and SDW channels 



In this subsection we derive the RG flows of density-wave channels, which have negligible overlap with the BCS 
channel. The channel is 

S DW = I d 2 xdt ^^^V tQ (ks;x) ^(ks + Qijas) 

S,T i 

x ^ (k T ; x) ^ (k T + Q 3 ; x) [8^ -8%]. (32) 

Here Qi,Qj are two of the four nesting vectors (Fig. |^). Due to the nesting, this channel flows logarithmically in 
scale length and causes a CDW or a SDW instability, depending on the sign of Fdw- The spin indices in Eq. (p3^) 
are appropriately symmetrized according to Fermi statistics.EJ Separating the interaction Fdw into decoupled CDW 
and SDW channels, 

1 f ' 

Sdw = 7j / d2x dt X] { F sdw Sq{S;x) ■ S Q (T;x) + F C dw n Q {S;x) n Q (T;x)} , (33) 

J S,T 
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where 



and 



S Q (S;x) = - V tQ (ks;x) o* ^(k s + Q;x) 



i Q (5;x)=^ta( ks;x ) </, Q (k s + Q;x) 



(34) 



(35) 



are the spin density and the charge density, respectively. The primed summation over S, T indicates that the sum is 
over the nested parts of the Fermi surface. First, we investigate the RG flow of the CDW term. For simplicity we 
drop the spin indices. From Eq. ([ll]) the bosonized form of the CDW action is 

1 ' r 

Scdw = »J2 / dt * x f cdw(S,T) ^(k s ;x) V(k s + Q,;x) ^(k T ;x) t/>{k T + Q r ,x) 

S,T J 

2 V (2tt) 2 



S,T ' 

AT^p) l^J dtdx F C dw(S,T) 



S,T ' 
47T, 



r v4tt 

x exp r~fT ^ ks + ; ^ + ^ kT + Q ; x ) - <^( k s ; 

where x denotes (x, t). The scaling of e c , x , t follows 

t — ► t'/s , xii — ► Xm /s , e c — > se^. , 



(k T ; x' 



>}■ 



(36) 



(37) 



where s is the scale of the length. Dimensional analysis shows that Scdw is scale invariant and thus marginal. 

Next the fields may be split into higher and lower-energy modes 4>(ks; x) — > 0'(ks; x) + ft(kg; x) where ft is the field 
of higher-energy modes whose momentum ranges X'(S)/2 < |p| < X(S)/2 where X'(S) = s' c /vf(S). The equal-time 
correlation function of h fields is 



/t /i v w, n « J, , ^ ns-x + 2 ?S /A(5) 

ft k s ; x ft k T ; = 5 SiT — In , . /ug x 

47r \ ns • x + 2i/X(b) 

= ; |n s x x| > 1 . 



|n s x x| < 1 



(38) 



We study how the partition function is renormalized after the fast momentum modes are integrated out: 

1 



-SCL 



(1-5, 



CDW 



c2 
, a CDW 



— ex P [ — S'cdw] 



(39) 



Here the average in the brackets is taken only over the high-energy fields ft and S' CDW is the renormalized CDW 
action. The third term on the right hand side (RHS) of Eq. ( p9|) gives the second order in the Fqdw correction to 
Scdw- Using the formula (e A e B ) = e (AB+(A 2 +B 2 )/2) and Eq. (|f), we obtain 



(Scdw) ~ (Scdw 



A 



(2^ 



S.T,U 



^2 dt d 2 x F 2 DW 



X(S} X(T) X 2 (U) 
16s 4 



r V 4tt 1 
exp j«-^-[<Xk s + Q t ;x) + 4>{k T + Q^x) - <f>{k s ;x) - <f>(k T ;x)]j 

J du d 2 y exp|^[(ft(ku;y,u) ft(k u; 0,u)) + (ft(ku + Q,; y, u) h{k v + Q ; ; 0, u))} } 



subleading terms 

4 / 



A 



(27T) 



Y, [dt d 2 x J2 f cdw ^pr x / du d(nu ' y) 

S,T J U 



4s 2 /X(U) 2 



[v F {U)u] 2 + (fiu -y) 



A(S) A(T) 



4s 2 

+subleading terms . 



exp |i— ^— [0(k s + Qi,x) + (j)(k T + Q^x) - </>(k s ;x) - ^(k T ;x)]j 



(40) 
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Here vf(U) denotes the Fermi velocity at the patch U. Replacing X(S) with sX'(S) [X'(S) is the new cutoff] and 
performing the integrals over ftu ■ y and u, we extract the correction to the CDW coupling function 



dFcDw(s) A 



E 3n f cdw ■ (41) 



dln(s) 2(2tt) 2 ^ v F {U) 

Therefore, we expect instability in this channel if Fqdw is strictly negative. This agrees with the result obtained by 
ShankarJa Fqdw < corresponds to a bare attractive interaction (U < 0) and the CDW instability leads the phase 
transition into a charge-density ordering with period (tt,tt) and (w,—n). 

Likewise, a similar RG transformation can be performed to obtain the flow equation for the SDW channel: 

1 



Ssdw =-- d'xdt }_^Fsdw Sq(S;x) • S Q (T;x) (42) 



S,T 



yielding 

dFsDw(s) _ A A 1 



d\n(s) (2tt)2 ^ v F (U) 



Fl DW (43) 



As expected there occurs an instability when Fsdw < 0, which corresponds to a bare repulsive interaction (U > 0). 
The SDW instability leads to an AF spin-density ordering with period (jr, 7r) and (tt, —n) since the on-site repulsive 
interaction induces an AF exchange interaction with the nearest neighbors. 

Since the RG equation for the BCS channel is of interest, we now discuss the overlap between the BCS and the 
SDW (CDW) channels. In Fig. | we show a typical scattering channel which belongs to both the BCS and the 
density-wave channel. In terms of the BCS, it can be written as 



Overlap = \ d2x dt ^( k S,k T ) ^"(ks^) ^ ("k S ; z) 



S,T ■ 

x^/3(-k T ;a;) ^ Q (k T ;^) <5 kT ,k s +Q • (44) 

Obviously there is a severe phase space restriction due to the factor of £k T ,k s +Q m the overlapped channel, which 
amounts to a factor of Aa -C 1. Therefore, the SDW (CDW) is practically decoupled from the BCS channel. The 
renormalization of the BCS channel comes exclusively from the other repulsive or AF exchange channels which we 
discuss in the next subsection. Here we can evaluate the energy scale of the SDW instability for future reference. 
From Eq. (|43]), we obtain the critical energy at which the SDW instability occurs: 



Esdw ~ E F exp 



1 



j50z In | \F SDW \ 



(45) 



B. Kohn-Luttinger effect 

The two classes of repulsive channels which renormalize the BCS are shown in Fig. |] and Fig. || and are given by 

S * = \ E E / * x dt F R ^ tQ (ks; x) ^(k s + Ap + Q i; x) 



S,T Ap 

7 7 

which is the zero momentum process (Fig. |j) and 



x^7 (kT . x ) ^ 5 (k T - Ap - Qi,x) [Sffi - , (46) 



S u = \Y.Y.J d2xdt F u ^ ta (ks; x) ^(k s + Ap + Q, t ; - 



S,T Ap 



x^(k T -x) M^t- Ap + Q^x) [8^-5%] , (47) 
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which is the 2kp — Q process allowed because 2Q is a commensurate wave-vector (Fig. ||). In bath channels we take 
the coupling functions to be momentum-independent and the spin indices are anti-symmetrized.Ej 

Now we expand the partition function to second order in Sr + Su, integrate out the high-energy modes, and 
study the resulting effective action generated by this procedure. In the course of this RG procedure, Sp and Su are 
renormalized as well as the BCS channel, and so the resulting flow equation will be a complicated set of coupled 
differential equations. We do not attempt to perform a complete analysis of the flow equations but we are only 
interested in the trend of initial growth of the BCS channel; the higher-order coupled equations are important only 
in the very-long-length scale. Thus we only study how Sr + Su renormalizes the BCS channel in the initial stages 
of growth near s = 1. The second-order correction to the BCS channel is calculated in the same manner as in the 
previous subsection: 



MSr + Su) 2 ) 



Ins (Fl, + 4F^)J2J2 I d 2 x dt^ a (k s ;x) ^(-k s ;x) 



S Ap ' 

X^kg + Ap + Q; x) ^ Q (-k s Ap Q; x) 

1 , /l + cos(e + |Ap|/2)\ 
in 1 1 



V2 

2(2tt) 2 cosAp/2 Vl - cos(e+ \Ap\/2)J 
In s 
~~2 



V2 



y\— 

4^ L(2tt) 2 



S.T 



In 





cos e\ 


1 1 




cos eJ 


+ d 




-kg; a;) 





non — BCS terms 



(48) 



The second term in Eq. (|48|) contributes to the s-wave BCS channel. Since it is repulsive even in the first order to 
begin with, the s-wave BCS channel is not a candidate for instability. The first term in Eq. (48) gives RG flow to all 
BCS symmetry channels since it has a strong wave-vector dependence. More succinctly, the RG equation of the BCS 
channel is given by 



^VBCs(ks,k T ) 
dins 



= {F 2 R + AF^) 
A 



V2 



1 



In 



(f 



cos(e+|Ap(S,T)|/2) 



2(2tt) 2 cosAp(5',T)/2 
^2 - ) rn V BC s(k s ,k u ) V BC s(k u ,k T ) 



cos(e+|Ap(S,T)|/2) 



(2tt)2 ^ VF (U) 



(49) 



if S and T are in the same quadrant on the Fermi surface. When S and T are in different quadrants, 
dVBcs(ks,k T ) A ^ 1 



dins 



(2tt) 



v F {U) 



Vbcs (ks.ku) Vbcs (ku,k T ) 



(50) 



Ap(S, T) is a vector along the nested Fermi surface which is shown in Fig. |[ For instance, Ap(S, T) = | (kg — kx) ■ x— 
7r| if ks and kx are both in the first quadrant on the {k x , k y ) plane. The expression of Ap(S l , T) in the other quadrants 
can be obtained by inverting the signs of ks or kx appropriately. 

Now we are ready to check which symmetry channels are attractive by evaluating the generalized Fourier coefficients 
of the double expansion shown in Sec. EI] . Among the basis functions we introduced in Eq. (p7|), we especially focus 
on the fundamental modes whose Fourier coefficients are presumably the most significant. After performing the 
generalized Fourier transform of the RHS of Eq. ( pl9| ) we find that there are two diagonal Fourier coefficients which 
flow to negative numbers among the fundamental modes. Their rate of change is evaluated as follows: 



dC 



dCf 
dins 



a 2 E 



1 



d Vbcs (k s ,k T ) 



S.T 



v F (k s ) up (kx) 



din { 



-0.2 (F 2 +4^)-|C 2 11 | 2 /(2 7 r) 2 , 

1 d Vbcs (kg,k T ) 



S.T 



v F {k s ) vp (kx) 



-0.7 (F 2 R + 4 Fy) 



|c oo,2 



din s 

2 



/(27T) 



a 2 (l,k s ) a 2 (l,k T ) 



6i(0,ks) 6i(0,k T ) 



(51) 



(52) 
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where C\ (s = 1) = C\ (s = 1) = initially. These are evaluated when e w 0.1. Here we find that the inequality 

,. 3 > M- holds for a sufficiently small e. From the Landau theorem we can infer that there can be a BCS 
instability in cither the Bi (d x 2_ J( 2-wave) or A 2 (eight-lobed s-wave) channel as long as the SDW instability does 

not set in. Furthermore, we can infer from the inequality hn^j > rfl 2 J that for a small enough e, d x 2_ y 2-wave 
superconductivity is preferred. The reason the inequality holds is that the d x 2_ y 2--wave basis function &i(0, k) has a 
larger overlap than 02(1, k) does with the density of states at the Fermi surface which peaks at the van Hove points. 

Of course, d x 2_ y 2-w&ve superconductivity is not the only possible ground state near half-filling. There is competition 
between the SDW and BCS instabilities since both are favored in the presence of a bare on-site repulsive interaction. 
Slightly below half-filling, the nesting is not perfect but nevertheless effective as long as the effective theory around 
the Fermi surface covers a significant part of the nested Fermi surface. As the energy cutoff around the Fermi surface 
is reduced via the RG towards e$, where e$ is the energy deviation of the actual Fermi surface from the nested one, the 
SDW and AF exchange channels sto». flowing whereas the BCS channel continues to be renormalized logarithmically 
as its flow is independent of nesting.El If initially the SDW coupling is small enough or if e$ is large enough so that 
it does not develop an instability until e c — + eg, or Esdw < £s, there never occurs a SDW instability. Although 
it is difficult to compute the critical energy, the BCS instability eventually sets in at a sufficiently low-energy scale. 
From the argument above, the superconducting order parameter presumably has d x 2_ y 2-w&\e symmetry. This result 
is consistent with the expectation that AF exchange channel is conducive to d K 2_ y 2-wav^superconductivity and 
qualitatively in good agreement with the parquet solution obtained for a flat Fermi surface.Ej 



V. SUMMARY AND CONCLUSIONS 



In this paper we investigated the BCS instability in a nested Fermion liquid with a bare repulsive interaction. 
When the ZS channel was incorporated we found that the low-energy fixed point was of a marginal Fermi liquid type. 
Although the Fermi liquid state broke down, multidimensional bosonization was still applicable. Large-momentum 
scattering channels were then included via the perturbative RG procedure in the boson basis. Among the logarith- 
mically renormalized channels, the SDW channel showed the quickest flow due to the nesting, but it had negligible 
overlap with the BCS channel. The other repulsive or AF channels gave rise to a renormalization of the BCS couplings. 
We found that the d x 2_ y 2-wa.ve channel was rendered the most attractive. Close to but less than half-filling, when 
no instabilities develop in the SDW channel (Esdw < £s), a d x 2_ y 2-wa,ve BCS instability sets in at sufficiently low 
temperatures. This is in agreement with the expectations of Monthoux, Balatsky, and PinesJ3 and the result obtained 
by Zhelenyak, Yakovenko, and Dzyaloshinskii.a 

The procedure presented in this paper is useful for extracting information about the low-energy fixed points of 
lattice systems. Once a realistic band structure of the system near the Fermi surface is given as the input, one 
can analyze the phase structure and the order parameter symmetry at low temperatures. As regards the cuprate 
superconductors, it would be interesting to study the RG of a three-band Hubbard model. For a detailed examination 
of the low-temperature phase structure of the materials, numerical works on the RG equation are needed. In this 
manner, one may be able to construct a phase diagram as a function of temperature and the filling factor (or doping). 
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FIG. 1. The low-energy effective theory in the vicinity of the Fermi surface. The thin momentum shell is coarse-grained into 
squat boxes of a height A and a width A. Q is the nesting vector. 

FIG. 2. A typical density- wave scattering channel. 

FIG. 3. A typical overlap between density-wave and BCS channel. 

FIG. 4. The repulsive or AF channels with a zero momentum transfer which logarithmically renormalize BCS channel. 

FIG. 5. The repulsive or AF channels with a momentum transfer 2Q which logarithmically renormalize BCS channel. 
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